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Suppression of high px hadron spectra in high-energy heavy-ion collisions at different energies 
is studied within a pQCD parton model incorporating medium induced parton energy loss. The 
Pt dependence of the nuclear modification factor Raa{pt) is found to depend on both the energy 
dependence of the parton energy loss and the power-law behavior of the initial jet spectra. The high 
Pt hadron suppression at y/s = 62.4 GeV and its centrality dependence are studied in detail. The 
overall values of the modification factor are found to provide strong constraints on the lifetime of 
the dense matter. 
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I. INTRODUCTION 

The discovery of jet quenching in central Au + Au col- 
lisions at the Relativistic Heavy-ion Collider (RHIC) at 
Brookhaven National Laboratory has provided clear ev- 
idence for the formation of strongly interacting dense 
matter. The observed jet quenching includes suppres- 
sion of single hadron spectra at high px [1,2], disappear- 
ance of back-to-back correlation of high px hadrons [3] 
and the azimuthal anisotropy of high px hadron spectra 
in non-central Au + Au collisions [4]. The absence of 
these jet quenching phenomena in d + Au collisions [5-8] 
proves that they are indeed due to final state interac- 
tion. Detailed analyses [9] further indicate that they are 
caused by parton energy loss as predicted [10-13]. Using 
the parton energy loss extracted from experimental data, 
one can conclude that the initial gluon density in central 
Au + Au collisions at \J~s = 200 GeV is about 30 times 
higher than that in a cold nucleus [14,15]. 

The observed strong suppression of high px hadrons 
at RHIC is in sharp contrast to the results of Pb + Pb 
collisions at the SPS energy. Data from the WA98 [17] 
experiment show no or little suppression of high px (up to 
about 4 GeV/c) pions [18,19]. Even if one takes into ac- 
count the possible uncertainty in the reference p + p data 
[20], the suppression allowed by the data is still signifi- 
cantly less than in Au + Au collisions at RHIC, while the 
total charged multiplicities or the inferred initial parton 
densities only differ by a factor of 2. This implies that 
additional physics is at play in the energy dependence of 
the suppression of single hadron spectra in high-energy 
heavy-ion collisions. It could be the energy dependence 
of the Cronin effect due to initial state multiple scat- 
tering, the thermalization time and finite lifetime of the 
dense matter that limits the parton energy loss. The en- 
ergy dependence of the hadron suppression at large px 
was predicted in Ref. [9] and also was recently studied 
in Ref. [16], with parton energy loss proportional to the 



observed hadron multiplicity. 

In this brief report, we explore the sensitivity of the 
final high px hadron suppression to the lifetime of the 
dense matter as well as the dependence on the collid- 
ing energy. We will study the high px hadron suppres- 
sion at yjl = 17.2, 62.4, 200 and 5500 GeV. We focus, 
in particular, on both neutral pion and charged hadron 
suppression at y/s = 62 .4 GeV and study the constraint 
on the lifetime by the measurement of hadron spectrum 
suppression. 



II. ENERGY DEPENDENCE OF HIGH P T 
HADRON SUPPRESSION 

We will use a LO pQCD model [18] to calculate the in- 
clusive high-pT hadron cross section in A + A collisions, 
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where a(ab — > cd) are elementary parton scattering 
cross sections and t A {b) is the nuclear thickness func- 
tion normalized to J d 2 bt A (b) — A. We will use a hard- 
sphere model of nuclear distribution in this paper. The 
K « 1.5 — 2 factor is used to account for higher order 
QCD corrections. The hadron is assumed to have the 
same rapidity as the parton, y = y c and its fractional 
momentum is defined as z c = px/ptc- The parton distri- 
butions per nucleon f a / A (x a , Q 2 , f) inside the nucleus are 
assumed to be factorizable into the parton distributions 
in a free nucleon given by the MRS D— ' parameteriza- 
tion [21] and the impact-parameter dependent nuclear 
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modification factor given by the new HIJING parame- 
terization [22]. The initial transverse momentum distri- 
bution gA^kx, Q 2 , b) is assumed to have a Gaussian form 
with a width that includes both an intrinsic part in a 
nucleon and nuclear broadening. Detailed description of 
this model and systematic comparisons with experimen- 
tal data can be found in Ref. [18]. 

The effect of parton energy loss is implemented 
through an effective modified fragmentation function 
[23], 



D h/c (z c ,Q 2 ,AE c ) = (l-e-^) 



-D° h/ M,Q 2 ) 



+(^) Z fD° h/g{z ' g ,Q 2 ) 



+ e-^Dl /c (z c ,Q 2 ). (2) 



This effective form is a good approximation to the ac- 
tual calculated medium modification in the multiple par- 
ton scattering formalism [24], given that the actual en- 
ergy loss should be about 1.6 times of the input value 
in the above formula. Here z' c — Pt/(pt c ~ AE C ), 
z' g = {AL/X)pt/AE c are the rescaled momentum frac- 
tions and AE C is the total energy loss during an average 
number of inelastic scatterings (AL/X). The FF's in free 
space D^ c (z c ,Q 2 ) are given by the BBK parameteriza- 
tion [25]. 

In this study, we assume a 1-dimensional expanding 
medium with a gluon density p g (r, r) whose initial distri- 
bution is proportional to the transverse profile of partic- 
ipant nucleons. The total energy loss for a parton prop- 
agating this medium is 



AE(b, r, 



dr — p g (T,b 7 r + nT), 



ToPo 



(3) 



according to recent theoretical studies [13,26,27], where 
po is the averaged initial gluon density at r in a central 
collision. (dE/dL}id is the average parton energy loss 
over a distance R A in a 1-d expanding medium with an 
initial uniform gluon density po. Similarly, the average 
number of scatterings along the path of parton propaga- 
tion is 



(AL/X) 



J Tn 



dTap g (T, b,r + nr) 



(4) 



Inclusion of transverse radial expansion generates a faster 
dilution of the gluon density relative to 1-d expansion, 
but also results in longer propagation time in medium. 
These effects offset each other and the final total energy 
loss in the cases of 1-d and 3-d expansion is found very 
similar [28] . With a hard-sphere nuclear distribution, the 
gluon density profile p g (r,r) can be expressed in a sim- 
ple form and the analytic expressions of AE(b, r, <f>) and 
(AL/X) are given in the Appendix. 

The energy dependence of the energy loss is parame- 
terized as 



= e Q (E/p - 1.6) 1 ' 2 /(7.5 + E/p) , (5) 

according to the numerical results in Ref. [29] in which 
thermal gluon absorption is also taken into account in the 
calculation of parton energy loss. Fit to the most cen- 
tral Au + Au collisions at y/s = 200 results in eo = 1-07 
GcV/fm, p = 1.5 GcV and A = l/ap = 0.3 fm. The 
corresponding energy loss in a static medium with par- 
ton density p over a distance Ra is [27] dE /dL = 
(R A /2T )(dE/dL) ld w 14 GcV/fm. This is about 30 
times higher than the parton energy loss in a cold nu- 
cleus [14]. 

At different energies, we assume that eo an d Po are pro- 
portional to the measured hadron multiplicity as given in 
Ref. [30] . Since there is no experimental measurement at 
the LHC energy yet, we will use the model calculation [22] 
to extrapolate and assume (N c h /dr]) 5500 /(N c h/dr]) 200 ~ 
2.4. 

To study the dependence of the high px hadron sup- 
pression on the colliding energy, let us assume first that 
the lifetime of the dense matter is larger than the system 
size. Shown In Fig. 1 are the the nuclear modification 
factors 



Rab(pt) 



do h AB 
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for charged hadron (solid lines) and neutral pions (dashed 
lines) in central Au + Au (Pb + Pb at the SPS energy) 
collisions at different energies, from SPS s/s = 17.2 GeV 
to LHC y/s = 5.5 TcV. Here, 

(^binary) = J d 2 bd 2 rt A (r)t A (\b - r\) (7) 

is the number of geometrical binary collisions at a given 
range of impact parameters. At the SPS energy, the ob- 
served nuclear modification factor in central Pb + Pb col- 
lisions is consistently about 1 due to strong Cronin effect 
via initial multiple parton scattering, leaving not much 
room for large parton energy loss [19]. We shall return 
to this point later. 

In central Au + Au collisions at RHIC, however, strong 
suppression of high pt hadrons is observed. This can be 
attributed to large parton energy loss that overcomes the 
modest Cronin enhancement as observed in d + Au col- 
lisions [5-8] and gives rise to the large hadron suppres- 
sion. The energy-dependence of the parton energy loss 
in Eq. (5) describes well the flat pr dependence of the 
nuclear modification factor R A a(pt) for neutral pions at 
large pr in Au + Au collisions at y/s = 200 GeV, as 
shown in the figure. For charged hadrons there are com- 
plications from other medium effects as we will discuss 
later. The flatness of the modification factor at this en- 
ergy is more clearly illustrated by the ratio of central to 
peripheral collisions [31,32]. Such a flat px dependence 
is actually a coincidence, as the combined effect of the 



2 



energy dependence of the parton energy loss in Eq. (5) 
and the power-law behavior of the initial jet spectra. Us- 
ing the same energy dependence of the parton energy 
loss but with a reduced amplitude due to smaller initial 
gluon density at yfs = 62.4 GeV, the nuclear modifica- 
tion factor is found to decrease with px and even becomes 
smaller than the modification factor at y/s — 200 GeV 
at high pt > 10 GeV/c. This is simply a consequence 
of the energy dependence of jet spectrum shape. The 
initial jet spectra at ^/s — 62.3 GeV are much steeper 
than those at 200 GeV. The same amount of energy loss 
leads to a larger suppression of the final hadron spec- 
tra at 62.3 GeV than at 200 GeV. As one increases the 
colliding energy, the power-law spectra for the initial jet 
production become flatter, and the same parton energy 
loss will lead to less suppression of the final hadrons. As 
shown in the same figure, the nuclear modification factor 
at the LHC energy y/s — 5.5 TeV is smaller than at 200 
GeV in the intermediate pt region, due to larger initial 
gluon density. However, the modification factor Raa{pt) 
increases with pj< due to the flatter power-law spectra of 
jet production at LHC. 




Pt (GeV/c) 

FIG. 1. Nuclear modification factors for charged hadrons 
(solid) and neutral pions (dashed) in 0-5% central Au + Au 
collisions at y/s = 17.3, 62.4, 200, 5500 GeV. The lifetime 
Tf — 0,1,2 fm/c (from top to bottom) is assumed for calcula- 
tions at y/s = 17.2 GeV. For other energies Tf is assumed to 
be larger than the system size. The STAR [31] and PHENIX 
[32] data are for central Au + Au collisions at y/s = 200 GeV. 

The parameterized energy dependence of the parton 
energy loss is in part due to the detailed balance effect 
in induced gluon radiation and absorption [27]. This ef- 
fect is most important in the intermediate pt region. In 
this region, one expects the jet fragmentation process to 
be modified by other non-perturbative processes such as 
parton recombination or coalescence [33-35]. The ob- 
served flavor dependence of the hadron suppression and 



of the azimuthal anisotropy clearly points to the effect 
of parton recombination that enhances both baryon and 
kaon spectra in the presence of dense medium. To include 
this effect in the current parton model, we have added a 
soft component to kaon and baryon FF's that is propor- 
tional to the pion FF with a weight ~ (Nbm(b, r))/[l + 
cxp(2pTc~ 15)]. The functional form and parameters are 
adjusted so that (K + p)/n ~ 2 at pt ~ 3 GeV/c in the 
most central Au + Au collisions at y/s = 200 GeV and 
approaches its p + p value at p T > 5 GeV/c. This gives 
rise to the splitting of the suppression factor for charged 
hadrons and tt° in the calculation. Because of the steeper 
power-law spectra of jet production at 62.4 GeV, the ef- 
fect of the non-perturbative parton recombination per- 
sists to higher px than in Au + Au collisions at 200 GeV. 
In this region of pr, the non-perturbative recombination 
effects dominate the nuclear modification of the charged 
hadron spectra. As a consequence, the px dependence 
of the modification factors Raa(pt) at y/s = 62.4 and 
200 GeV are similar. They only diverge at high pt > 8 
GeV/c where the recombination effects are negligible at 
both energies. Since the soft components due to parton 
recombination are not closely related to medium ther- 
malization, they will still contribute to the final hadron 
spectra in peripheral A + A and p+A collisions. This will 
lead to flavor dependence of the Cronin effect in p + A 
collisions [36]. 

III. EFFECT OF FINITE LIFETIME 

In the calculation of the parton energy loss in Eqs. (3) 
and (4), the upper limit of the path integral should be 

Tmax = min(AL,r/), (8) 

where r/ is the lifetime of the dense matter before 
breakup, AL(b, f, 4>) is the distance the parton, produced 
at r, has to travel along n at an azimuthal angle <f> relative 
to the reaction plane in a collision with impact-parameter 
b. According to the analyses of experimental data on high 
Pt hadron suppression and suppression of away-side cor- 
relation in Au + Au collisions at y/s = 200 GeV, the 
extracted energy loss points to an initial gluon density of 
about 30/fm 3 at an initial time r = 0.2 fm. Given the 
measured transverse energy per charged hadron of 0.8 
GeV [37], this gives a lower bound on the initial energy 
density of about 25 GeV/fm 3 . In 1-d expansion with 
the equation of state of an ideal fluid, the energy den- 
sity decreases with time, e(r) = e^Ti/r) 4 / 3 . Assuming 
the 1-d hydrodynamics expansion starts at Ti = 1 fm/c 
(free-streaming before that), the lifetime of the plasma or 
the duration for the parton energy loss should be about 
Tf ~ 5 fm, before the phase transition with a critical 
energy density e c ~ 1 GeV/fm 3 . The early stage of the 
mixed phase or crossover could also contribute to the jet 
quenching and thus extend the effective time duration 
for parton energy loss. When this time is larger than the 
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average path length, the total parton energy loss is then 
limited only by the system size. In the previous analysis 
[14,29], such an assumption for central Au + Au collisions 
is justified given the high initial energy density. 

In central Pb + Pb collisions at the SPS energy y/s = 
17.3, the rapidity density of charged particles is about 
half of that in central Au+Au collisions at y/s — 200 GeV 
[30] . The average transverse energy per charged particle 
is about the same at SPS and RHIC energy [37]. One can 
then assume the initial energy density at SPS to be half 
of that in central Au + Au collisions at y/s = 200 GeV. 
The lifetime of the plasma in central Pb + Pb collisions 
at SPS, t s ~ 2 — 3, if it were formed, is then consider- 
able shorter than that at the highest energy of RHIC. 
The longer thermalization time at SPS could make the 
effective lifetime even shorter. To demonstrate the de- 
pendence on the lifetime, we show in Fig. 1 the nuclear 
modification factor of tt° at the SPS energy with short 
lifetime r/ = 0,1,2 fm/c. The modification factor is 
quite sensitive to the lifetime Tf. Even if one takes into 
account the uncertainty in the reference p + p data [20] , 
which gives a systematic error of about a factor of 2, the 
data still point to a short lived dense system with rt < 2 
fm/c. This might also explain the observed elliptic flow 
i>2 at SPS that is much smaller than the hydrodynamic 
limit [38] . On the other hand, the longer lifetime of the 
plasma allows development of a full hydrodynamic flow, 
giving rise to a large V2 that saturates the hydrodynamic 
limit [39]. 
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FIG. 2. The centrality dependence of nuclear modification 
factors for charged hadrons (solid) and neutral pions (dashed) 
in Au + Au collisions at yfs — 62.4 GeV. The lifetime of the 
dense medium is r/ = 10,4,3,2,1,0 fm/c (from bottom to 
top). For central collisions results with only large values of 
Tf are presented. The dot-dashed line in the first panel is for 
it in 0-10% central collisions with Tf = 10 fm/c 



To study the sensitivity of the high pt hadron suppres- 
sion to the lifetime of the plasma at y/s = 62.4 GeV, we 
show in Fig. 2 the nuclear modification factors for both 
charged hadrons and neutral pions with different values 
of Tf in Au+Au collisions at y/s = 62 A GeV. The hadron 
suppression in the large px region in the most central col- 
lisions is very sensitive to the lifetime of the plasma in 
this calculation. In peripheral collisions, the small size 
of the dense medium limits the parton energy loss. As a 
result, the hadron suppression is only sensitive to values 
of Tf that are smaller than the average medium size. In 
reality, the values of r/ should decrease from central to 
peripheral collisions. 

The recent experimental results from PHOBOS [40] 
on nuclear modification factors for charged hadrons in 
Au + Au collisions at 62.4 GeV only extend to px ~ 
4 GeV/c. In this region, the suppression of charged 
hadron is indeed much smaller than at 200 GeV. How- 
ever, charged hadrons in this region are also dominated 
by non-perturbative recombination effects, though our 
results are still sensitive to the lifetime. Experimental 
measurements of ir° and high px charged hadrons, both 
are less influenced by the parton recombination effect, 
should provide more stringent constraints on the lifetime 
of the dense matter. 



IV. SUMMARY AND DISCUSSION 

Within a parton model incorporating medium induced 
parton energy loss, we have studied in this brief report 
the suppression of inclusive hadron spectra at high px 
in heavy-ion collisions at different energies. We found 
that the px dependence of the nuclear modification fac- 
tor Raa(j>t) is determined by the energy dependence of 
the parton energy loss and the power-law behavior of the 
initial jet spectra. With the onset of parton energy loss 
and the change of the power-law jet spectra, the px de- 
pendence of the modification factor changes from mono- 
tonic decrease at y/s = 62.4 GeV to monotonic increase 
with px at the LHC energy yfs — 5.5 TeV. The flat px 
dependence observed at y/s = 200 GeV is just a coinci- 
dence. 

We also studied the sensitivity of the hadron suppres- 
sion factor to the lifetime of the plasma or the duration 
of parton energy loss. We found that the hadron sup- 
pression factor at intermediate and large px is sensitive 
to the lifetime if it is comparable or smaller than the sys- 
tem size. The experimental measurements could provide 
important constraints. Together with the measurement 
of elliptic flow, which is also sensitive to the lifetime and 
thermalization time of the plasma, one can gain addi- 
tional information on the dynamic evolution of the pro- 
duced quark-gluon plasma. 

One can also calculate the back-to-back dihadron cor- 
relation at different energies. We find that suppression 
of the back-to-back dihadron correlation at 62.4 GeV in 
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central Au + Au collisions is almost identical to that at 
200 GcV. This is partly due to the trigger bias that se- 
lects dihadron production close to the surface and results 
in completely suppression the back-side jets that traverse 
the whole length of the dense matter. The suppression 
due to kr broadening of initial multiple parton scattering 
is also independent of the colliding energy. 
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V. APPENDIX 

In this appendix, we give the basic analytic formula for 
calculating the path integral in the parton energy loss in 
Eqs. (3) and (4), assuming a hard-sphere nuclear distri- 
bution. 

Given two overlapping nuclei as illustrated in Fig. 3, 
we want to calculate a path integral over the path AL. 
Let n and r 2 be the radial coordinates of the jet produc- 
tion point as measured from the center of the two nuclei. 
For given b and r\ 



r 2 = \jb 2 + r 2 - 2bn cos<£ 6 (9) 
For a hard-sphere distribution, r x < Ra and r 2 < Rb- 




FIG. 3. Transverse geometry of two overlapping nuclei 

The jet travels at an angle (j> p with respect to the 
impact-parameter b (which determines the reaction plane 
with the beam direction). <pb is the angle between b 
and f\. The distance the jet travels is then the AL = 
min(ri,r 2 ); 



T\ = \J R\ — r\ sin 2 <jV — r\ cos <j) r ; 

r 2 = \J R 2 B — r 2 + (ri cos <j) r — b cos 4> p ) 2 



(ri cos ( 



6cos( 



(10) 



(11) 



where </> r = cf) p — <t>b- These are solutions of 

Ra = \n +nn|, (12) 
Rb = \r 2 - b + nr 2 \ (13) 
If b< \R A -R B \ 

AL = \l R2 min - r i sin2 06 - n COS <j> b , (14) 

where R m { n — Miu(Ra, Rb)- We also define i? max = 
Ma,x(R A ,RB)- 

Assuming that the soft gluon density is proportional 
to the number of participant nuclcons, it is then given by 



p g (T,b,f) = 



tqPq 



t 2c A bRu 



R" 



t A (r)6(R B -\b- r\) 



B 



■t B (\b-f\)e(R A -r) 



(15) 



where c AB = 1 - (1/2) (1 - R^JR^f' 2 and p is de- 
fined as the averaged gluon density in central collisions 
(b = 0) at an initial time r : 



Po 



irR 



2 

min 



d 2 rp(T Q ,r,b = 0). 



^) = ^l-r 2 /R% 



Using the nuclear thickness function defined as 

3 A_ 

l 2 

l A 

the gluon density is then 
3t Po 



(16) 



(17) 



p g (r,b,r) = 



AtcabR 11 



6{R A -r)^R 2 B -\b-r\ 2 



9(R B -\b-r\)jR A 



(18) 



According to Eq. (3), the total energy loss along the 
path AL is 



dE 

AE(b,r,4>) = {—) ld 
dE 



AL+r 



dr- —P 9 (t 7 b,r + nr) 



dL' lu 4cA B R 
^R\-(r + nr) 2 



L 



TOPO 
t +AL 



dr 



T - T Q 



+ ^R 2 B -\b-{r + nr)\ 2 



The average number of scatterings is 
(AL/X) = / dTap g (T,b,f+nT) 

J Tn 



(19) 



'T 

3 t 



f r„+AL ^ 
4A i? min Jr T 

+ \J R 2 B - \b- (r + nr)| 2 



R\ - (r + nT) 2 

(20) 
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The above integrals can be completed analytically. 
The following arc some basic integrals: 



/ 



dr^R 2 -(r + nr) 2 = T ^ ' — \J R 2 — (r + nr) 2 

R 2 — r 2 + (r ■ n) 2 r + r '■ n 

H arcsm 



^R 2 -r 2 + (r-n) 2 



(21) 



J ^-^R 2 -(f+nr) 2 = ^R 2 -{r + n T ) 2 



- (r • n) arcsin 



r ■ n + t 



y/R 2 -r 2 + {r-nf 



-sjR 2 -r 2 {log [R 2 -r 2 - (f ■ n)r 

- \/R 2 -r 2 ^jR 2 - (f+nr) 2 -logrj (22) 
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